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ABSTRACT 
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We obtained relatively high signal-to-noise X-ray spectral data of the early massive star t Sco (BO. 2 V) 
with the Suzaku XIS instrument. This source displays several unusual features that motivated our study: 

(a) redshifted absorption in UV P Cygni lines to approximately +250 km s _1 suggestive of infalling gas, 

(b) unusually hard X-ray emission requiring hot plasma at temperatures in excess of 10 MK whereas 
most massive stars show relatively soft X-rays at a few MK, and (c) a complex photospheric magnetic 
field of open and closed field lines. In an attempt to understand the hard component better, X-ray data 
were obtained at six roughly equally spaced phases within the same epoch of r Sco's 41 day rotation 
period. The XIS instrument has three operable detectors: XIS1 is back illuminated with sensitivity down 
to 0.2 keV; XISO and XIS2 are front illuminated with sensivitity only down to 0.4 keV and have overall 
less effective area than XIS1. The XISO and XIS3 detectors show relatively little variability. In contrast, 
there is a * 4<x detection of a « 4% drop in the count rate of the XIS 1 detector at one rotational phase. In 
addition, all three detectors show a*3% increase in count rate at the same phase. The most optimistic 
prediction of X-ray variability allows for a 40% change in the count rate, particularly near phases where 
we have pointings. Observed modulations in the X-ray light curve on the rotation cycle is an order of 
magnitude smaller than this, which places new stringent constraints on future modeling of this interesting 
magnetic massive star. 

Subject headings: stars: r Sco; stars: early-type; stars: magnetic fields; stars: winds, outflows; X-rays: stars 



1. Introduction 

The B0.5V star r Sco has a storied history in the 
study of massive stars. A northern hemisphere object 
of early spectral type that is relatively bright and close 
has made r Sco a favored target for multi-wavelength 
studies. It was a primary target in the early days 
of ultraviolet astronomy (e.g., Rogerson & Lamers 
1976; Rogerson & Upson 1977) and was found to 
be a relatively hard X-ray source as compared to its 
peers (MacFarlane & Cassinelli 1989; Cassinelli et al. 
1994). It also displays redshifted absorption in UV 
P Cygni lines of Ovi and Nv (Lamers & Rogerson 
1978). These features combined with growing evi- 
dence of clumped wind flows motivated Howk et al. 
(2000) to develop a model of infalling wind clumps 
to explain the observed redshifted absorptions and as a 
way of producing high velocity shocked gas to account 
for significant amounts of high temperature gas in ex- 
cess of 10 MK in the system. X-ray data from Chandra 
(Cohen et al. 2003) and XMM-Newton (Mewe et al. 
2003) confirmed the unusual hardness of this source, 
revealing also that the emission lines are relatively nar- 
row. Since then, advances in observational spectropo- 
larimetry have allowed for more detections and more 
stringent constraints of magnetism in massive stars 
(e.g., Donati et al. 2006a; Hubrig et al. 2008; Alecian 
et al. 2008). Donati et al. (2006b) reported on the dis- 
covery of a surprisingly complex photospheric mag- 
netic field in r Sco, a topology that sports sectors of 
open field lines and confining loops that deviate sig- 
nificantly from a simple dipole morphology. 



All of these observational features together have yet 
to be explained by a single unifying model. The emis- 
sion lines resolved by Chandra span a large range of 
energies, and yet all of the lines are quite narrow, with 
typical HWHM velocities of 300 km s . From Gaus- 
sian fits, the peak emissions for an ensemble of lines 
show little net velocity shift (Cohen et al. 2003). Both 
of these features suggest minimal photoabsorption by 
the wind of t Sco (e.g., Owocki & Cohen 2001). It 
may also be that the wind is relatively porous. (Dis- 
cussions of the impact of clumping and porosity for 
X-ray emissions and line profile shapes can be found 
in Feldmeier, Oskinova, & Hamann (2003), Oskinova, 
Feldmeier, & Hamann (2004), and Owocki & Cohen 
(2006). 

An understanding of the X-ray hardness and line 
profile properties is critical for producing a quantita- 
tive physical model for t Sco's wind structure. One 
way to further constrain the modeling effort is to ana- 
lyze variability of t Sco. Observing time of t Sco us- 
ing Suzaku was awarded to obtain multiple pointings 
with rotational phase to constrain the relative contri- 
butions of confined plasma to wind shocked plasma in 
terms of the X-ray emission measure and plasma tem- 
perature. Our experiment was designed to take advan- 
tage of the following facts. First, r Sco is a slow ro- 
tator with a rotation period of approximately 41 days. 
It is also relatively bright in X-rays, and so exposures 
for the different pointings are quite short as compared 
to the rotation period. In addition, the star is viewed 
nearly edge-on with respect to the rotational equator 
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such that the closed magnetic loops, which are asym- 
metrically distributed around the star, are occulted dur- 
ing a rotation cycle. In fact, Donati et al. (2006b) pre- 
dicted that the X-ray emissions of t Sco could vary by 
as much as 40% with rotational phase. 

Given that the field topology has remained stable 
over the last decade (Donati & Landstreet 2009), mul- 
tiple pointings - even if at different epochs - have the 
potential of determining hot plasma components in the 
closed magnetic loops as compared to the larger scale 
wind through an analysis of X-ray variability. This is 
a technique that has been used successfully in another 
magnetic massive star, 6 l Ori C. Gagne et al. (1997, 
2005) reported on multiphase X-ray spectral data from 
this young 05.5V star. In that case the surface field 
has kiloGauss values (Donati et al. 2002). The field 
is sufficiently strong in relation to its wind kinetic en- 
ergy density that substantial confinement of the wind 
flow in the form of an oblique magnetic rotator results, 
leading to significant X-ray variability with rotational 
phase (Babel & Montmerle 1997a). Gagne et al. re- 
ported changes in the soft and hard component emis- 
sion measures (EM) at about 25%. 

The method has also been used successfully for 
constraining emission models in the case of the mag- 
netic early B star /3 Cep (Favata et al. 2009). This is a 
Bl Hie star that is slowly rotating (vsim * 20 km s -1 , 
Abt et al. 2000). Despite a relatively strong surface 
field of about 360 G (Henrichs et al. 2000), a Chandra 
and XMM-Newton study by Favata et al. (2009) with 
pointings at four different rotational phases of /3 Cep 
failed to detect the expected modulation of X-ray emis- 
sion. The absence of significant variability strongly 
constrains the distribution of magnetically confined 
hot plasma in that system. 

In the case of t Sco, the surface field is slightly 
stronger at 500 G than for /3 Cep and considerably 
weaker than for 6 1 Ori C. Cohen et al.'s (2003) Chan- 
dra study of t Sco claimed variability at the 10% level 
in the O vm 18. 97 A line (0.654 keV) from an effective 
exposure of 72 ksec that was distributed over an in- 
terval of 120 ksec, or 3% of the rotation period of the 
star (see Fig. 4 of that paper), r Sco was also a regular 
target for the International Ultraviolet Explorer (IUE), 
and the P Cygni lines of Nv, Siiv, and Civ show vari- 
able structure with rotational phase. So, there is moti- 
vation for seeking multiple pointings toward t Sco in 
an attempt to exploit the "screening" effects afforded 
by stellar occultation in hopes of ascertaining the lo- 
cation in radius and azimuth of the hot plasma com- 



ponents in this magnetized massive star. A description 
of the observations and data reduction is presented in 
Section 2, followed by an analysis of the multiphase 
spectral data in Section 3. A discussion of our results 
is given in Section 4 with concluding remarks given in 
Section 5. 

2. Observations and Data Reduction 

The joint Japan/US X-ray astronomy satellite 
Suzaku observed r Sco in August and September 2008 
at six roughly equally spaced occasions. The pointings 
with 12-15 ksec exposure time were spaced approxi- 
mately 7 days apart to sample one full rotational cycle 
of the star which has P rot = 41. 033 days (Donati et al. 
2006b). Table[T]lists the obsids and dates of our point- 
ings along with the exposure times in columns 1, 2, 
and 4. Other columns in this table will be described 
later. 

Suzaku carries four X-ray Imaging Spectrometers 
("XIS", Koyama et al. 2007). The field-of-view (FOV) 
for the XIS detectors is 17' x 17'. One of the XIS 
detectors (XIS1) is back-side illuminated (BI) and the 
other three (XIS0, XIS2, and XIS3) are front-side il- 
luminated (FI). The bandpasses are ~0.4- 12keV for 
the FI detectors and ~ 0.2- 12keV for the BI detector. 
The BI CCD has higher effective area at low energies; 
however its background level across the entire band- 
pass is higher as compared to the FI CCDs. The XIS2 
detector was not used in our t Sco observations ow- 
ing to technical difficulties. Cleaned event files were 
used to extract the spectra for the science analysis. Re- 
sponse matrices were generated for each detector and 
each pointing using the latest versions of software and 
the most recent calibration files. 

The angular resolution of the X-ray telescope on- 
board Suzaku is « 2'; and so in the XIS imaging, 
t Sco was not distinctly resolved from the nearby X- 
ray source 1RXS J 163556.0-281 149. According to 
the recommendations given in the Suzaku ABC guid^H 
for spectral analysis, the spectra were extracted from 
a 260 arcsec region. Inspection of XMM-Newton im- 
ages shows that three serendipitous X-ray sources are 
present in the spectrum extraction region. The count- 
rates of these sources are a few hundreds of times 
lower than the count rate of r Sco. No variations in 
the X-ray emissions of the two sources located within 
2 arcmin of t Sco were reported from previous XMM- 



http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
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Table 1 

SUZAKU Observations of t Sco 



ObsID 


HJD-2454696 
(days) 


Phase 

4> 


Exposure 
(ksec) 


Ct 


XIS0 

Cs 


Ch 


Ct 


XIS1 

C s 


C H 


C T 


XIS3 

Cs 


C H 


403034050 
403034060 
403034010 
403034020 
403034040 
403034030 


0.31458 
7.19816 
14.04367 
21.07501 
27.42073 
33.41794 


0.168 
0.336 
0.503 
0.674 
0.829 
0.975 


16.3 
15.0 
14.6 
14.5 
14.0 
12.3 


0.806 ± 0.008 
0.808 ± 0.008 
0.825 ± 0.008 
0.829 ± 0.008 
0.810 ±0.009 
0.849 ±0.011 


0.306 ± 0.005 
0.317 ± 0.005 
0.329 ± 0.005 
0.325 ± 0.005 
0.326 ± 0.006 
0.338 ± 0.006 


0.500 ± 0.006 
0.491 ±0.006 
0.496 ± 0.006 
0.504 ± 0.006 
0.484 ± 0.007 
0.511 ±0.008 


1.554 ±0.010 
1.556 ±0.011 
1.498 ±0.011 
1.585 ±0.011 
1.562 ±0.012 
1.608 ±0.015 


0.859 ± 0.008 
0.862 ± 0.008 
0.820 ± 0.008 
0.875 ± 0.008 
0.862 ± 0.009 
0.890 ±0.011 


0.695 ± 0.007 
0.694 ± 0.007 
0.678 ± 0.008 
0.710 ±0.008 
0.700 ± 0.008 
0.719 ±0.010 


0.827 ± 0.008 
0.841 ± 0.008 
0.834 ± 0.008 
0.829 ± 0.008 
0.820 ± 0.009 
0.849 ±0.010 


0.324 ± 0.005 
0.334 ± 0.005 
0.339 ± 0.005 
0.333 ± 0.005 
0.335 ± 0.006 
0.345 ± 0.007 


0.503 ± 0.006 
0.507 ± 0.006 
0.494 ± 0.006 
0.496 ± 0.006 
0.485 ± 0.007 
0.504 ± 0.008 






Energy (keV) 

Fig. 1. — The spectrum of r Sco obtained by merg- 
ing the XIS1 spectra obtained on dates HJD 2454696, 
2454703, and 2454710 (see Tab.Q}, along with a best 
fit model (see text). The vertical error bars correspond 
to 10cr. Although no spectral lines are resolved, a 
number of typically prominent emission lines are iden- 
tified across the top of the figure. 

Newton or Chandra observations (Mewe et al. 2003; 
Cohen et al. 2003). Therefore, we believe that the 
extracted spectra are truly representative of the X-ray 
emission of r Sco. However the intrinsic variability of 
the serendipitous sources cannot be excluded a priori. 

The background was extracted from an annulus 
around the source for all six data sets. The presence 
of faint point sources in the background regions can- 
not be excluded (as seen in the XMM-Newton images). 
However, checks were performed using various back- 
ground regions selected in different parts of the Suzaku 
FoV, and it was determined that the influence of the 
background on the r Sco spectrum is negligible. 

3. Analysis 

3.1. X-ray spectra 

Spectra that were obtained at different phases were 
found to be quite similar. The spectral fitting software 
xspec (Arnaud 1996) was used for modeling the data. 
Based on chi-square fits, there is little evidence for sta- 
tistically significant differences between the six spec- 
tra obtained at the different pointings. Given the simi- 
larity of the spectra, we merged the first three spectral 
data sets obtained with XIS 1 to achieve higher signal- 
to-noise. The resulting spectrum is shown in Figure Q] 



The model spectral fit for this merged set is quite sim- 
ilar to fits obtained for the spectra from the individual 
pointing. 

A good fit to the spectrum was obtained with a four 
temperature thermal plasma model (vapec) corrected 
for interstellar absorption (tbabs). Abundances were 
taken as solar (Asplund et al. 2009), except for C at 
0.5 solar, N at 1.4 solar, and oxygen at 0.7 solar. These 
values are similar to the abundances derived by Morel, 
Hubrig, & Briquet (2008) in their study of UV and op- 
tical spectra of r Sco. Also, a ratio N/O of twice so- 
lar is in good agreement with the analysis of XMM- 
Newton spectra of tSco by Mewe et al. (2003). In 
fact, the good agreement between Suzaku and XMM- 
Newton spectra show that the Suzaku data are not seri- 
ously affected by the unresolved X-ray sources in the 
vicinity of t Sco. These assumptions lead to fits that 
have good statistics and provide a convenient frame- 
work for characterizing and comparing spectra at dif- 
ferent rotational phases. 

The best fit parameters found from fitting the spec- 
tra are shown in Table[2] The hydrogen column density 
was fixed at a value of A^h = 3 x 10 20 cirT 2 from Diplas 
& Savage (1994). It is interesting to compare the tem- 
perature components from our Suzaku fits to those of 
previous studies. An ASCA spectrum reported by Co- 
hen et al. (1997) was fit with a three temperature com- 
ponent model. The EM-weighted average temperature 
of that fit was (kT) asca ~ 1 keV. Ours is much lower at 
(kT) Suz * 0.5 keV. The XMM-Newton study by Mewe 
et al. (2003) is (kT} XM M ~ 0.7 keV from their RGS 
data. However, those authors find evidence for a hotter 
component from the EPIC spectrum. In their Chandra 
study of t Sco, Cohen et al. (2003) do not present tem- 
perature component fits; however, they do comment 
that there is little evidence for hot plasma in excess 
of 27 MK that was seen by ASCA. Following, they 
note that ASCA has greater sensitivity at higher ener- 
gies than does Chandra. In their study of differential 
emission measures for hot massive stars, Wojdowski & 
Schulz (2005) note that the peak of the EM distribution 
for t Sco is around 10 MK, or 0.8 keV. 

The XIS detectors have even greater area than 
did the ASCA-SIS detector, and yet the XIS1 spec- 
trum fails to show a significant contribution from the 
Sxv emission line. The total EMs derived from the 
ASCA data and from our Suzaku data are actually 
quite comparable, with values of « 49 x 10 53 cm -3 
and a 40 x 10 53 cm 4 , and yet it appears that the 
characteristic temperature during our pointings with 
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Suzaku is about half as large as when the ASCA data 
were obtained. 

3.2. Broadband variability 

Our primary objective was to search for evidence 
of stellar occultation effects on the X-ray emissions in 
an attempt to constrain the location where the soft and 
hard components are formed. Since the spectra show 
no obvious variations between pointings, we have pro- 
duced total count rate light curves. With a relatively 
complex distribution of closed and open field lines, we 
expected to observe variability with rotational phase. 

An ephemeris for the rotation of the field configu- 
ration with Heliocentric Julian date (HJD) is given by 
Donati et al. (2006b) as 



T(HJD) = T + PE, 



(1) 



where r (HJD) = 2453193.0, the period P = 41.033 + 
0.002 days, and E is the cycle value. All of our data 
fell within the same epoch of rotational cycle 1255. 
Rotational phase values within this cycle number are 
listed in column 3 of Table Q] Note that Donati & 
Landstreet (2009) show more recent magnetic maps 
for r Sco (see their Fig. 4). The most recent of those 
was taken at about 50 days before the Suzaku point- 
ings. The field has remained stable since the discov- 
ery report of Donati et al. (2006b), indicating that any 
accumulated systematic errors in rotational phases for 
the Suzaku data are less than 2%. 

Spectra for the three detectors with each pointing 
were summed to obtain total bandpass count rates, CV, 
in counts per second. These are listed in the last three 
columns of Table Q] and plotted in Figure [2] The hor- 
izontal lines in each panel represent the mean of the 
data, as computed via 



C, 



(2) 



where {en} are the error values given in Table Q] Note 
that the error in the mean is given by 



(3) 



for six data points obtained with each detector. The 
resultant means and their errors are given in column 2 
ofTable[3] 

The hashed regions in Figure [2] represent the rms 
errors cr^ for each data set as given by 



0.B6 
0.84 
0.82 

8 

(/) 
a, 
o 
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0.84 
0.82 
0.8 
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Fig. 2. — Total count rates from TableQjplotted against 
rotational phase. The colored hashed region represents 
the rms error of the data, as described in the text. 



(4) 



Given that the errors in the data are all very similar, the 
rms error closely represents the values of the individ- 
ual measures that determine the respective means. 

Although OB stars do show hot gas components, 
the hard emission is not seen with the relative strength 
at which it appears in r Sco (exceptions being OB stars 
that are strongly magnetic or in colliding wind binary 
systems). Thus, it is natural to associate the hot com- 
ponents seen from t Sco with processes involving the 
star's magnetism. Then the dominant soft emission 
most likely arises from nominal wind shock phenom- 
ena as occurs in other OB winds. Since an analysis of 
He-like ions that form at high temperatures indicate a 
formation location of about lR t above the photosphere 
(Cohen et al. 2003), we therefore anticipated that pe- 
riodic rotational modulation would affect the hard X- 
rays owing to stellar occultation, whereas any variabil- 
ity of soft X-rays would be mostly stochastic (e.g., Os- 
kinova et al. 2001). Consequently, we have defined 
"soft" and "hard" bands as illustrated in Figure [3] The 
upper panel displays the area sensitivity of the BI de- 
tector (XIS 1 ) as the upper curve in red and of the FI 



6 



Table 2 
Spectral Fit Parameters 



Component 


kT 


EM 




(keV) 


(10 53 cm- 3 ) 


#1 


0.11 +0.01 


17.0 + 2.6 


#2 


0.34 + 0.01 


10.4 + 0.5 


#3 


0.71+0.01 


7.2 + 0.3 


#4 


1.52 + 0.06 


5.2 + 0.3 



detector as the lower curve in blue. The middle panel 
shows the percent difference of these. For reference 
the lower panel displays a sample of the XIS1 spec- 
tra of t Sco (brighter points with black dots) and XIS3 
spectra (fainter with green dots). The vertical dashed 
line demarcates the soft and hard bands. Its placement 
at 0.95 keV was motivated by the following, (a) The 
boundary falls roughly midway along the decline in 
relative sensitivity, from the FI chip being about 75% 
more sensitive than the BI chips in the soft band to only 
25% on average in the hard one. (b) There are roughly 
equal counts above and below the boundary, so that 
count rate errors will be similar. And (c) the bound- 
ary falls at energies just below several strong emission 
lines. 

With these two bands, we introduce the hardness ra- 
tio HR = Ch/Cs as a ratio of the hard band count rate 
to that of the soft band. Figures 0}j6] plot HR values 
against the total count rate for each detector. An error 
ellipse is shown for reference. This is the rms error el- 
lipse with a center that is the mean of the HR and total 
count rate. Vertical and horizontal lines that are tan- 
gent to the error ellipse are shown as a visual reference. 
A light blue line connects the points in chronological 
and phase order, with the first (#1) and last (#6) point- 
ings indicated, except for XIS1 where several points 
fall too close together to label clearly which is first. 

Table [3] summarizes the mean values C av and HR dv 
derived from the data and the associated errors in the 
means. Also tabulated are reduced chi-square values 
xl with v = 5 for the null hypothesis of no intrinsic 
source variability. Recall that there is a 7.5% chance 
of achieving xl — 2 in the null hypothesis, but only 
a 1% probability of achieving ;f; = 3. It thus seems 
likely that we have detected variability of the total 
count rate in the XIS 1 detector, but less so for XIS0 
or XIS3 (however, see the discussion below about the 



high count rates for the last phase pointing). Despite 
variability in the count rate, the hardness ratios for 
XIS1 are remarkably stable. The spread in HR values 
for XIS0 and XIS3 is larger than for XIS1, although 
this greater spread is not very statistically significant. 

A primary culprit for the variability detected in 
XIS 1 is the third pointing near phase (f> « 0.5 as seen 
in Figure |2] Here C shows a 3.7% drop, whereas both 
XIS0 and XIS3 give no indication of being much dif- 
ferent from the mean. This could indicate variability 
at the softer energies detectable by XIS1, where XIS0 
and XIS3 have less areal sensitivity, especially below 
0.5 keV (see Fig. [3). The suggestion is made more 
plausible by recognizing that in Figure [5] the lowest 
count rate pointing of XIS 1 corresponds to the highest 
HR value. Since HR is a ratio of hard to soft, a reduc- 
tion in the soft count rate would increase the value of 
HR. 

We have made a comparison of the XIS 1 spectrum 
for the low count rate point at <p * 0.5 against spec- 
tra obtained at other phases. By binning XIS1 spectra 
for different pointings and overplotting them against 
the low count rate pointing at = 0.50, a systematic 
deficit in the count rate between about 0.8-0.9 keV is 
revealed. Thus the reduced count rate for <p = 0.50 
along with a slight increase in HR is not simply a 
"drop-out" in emission within a narrow energy band, 
but reflects a somewhat distributed reduction in emis- 
sion that lies below 1 keV near the peak of the spec- 
trum. Note that we also considered variability in the 
EMs of the four temperature components. In doing so, 
the four values of kT were held fixed, as well as the 
hydrogen column density. Indeed, the EMs of the two 
lower temperature components are more variable than 
those of the two high temperature ones, as expected. 

Interestlying, the ASCA data that shows evidence 
of a considerably higher temperature component than 
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Table 3 

Statistical Analysis of the X-ray Data 



Detector 


r 

(cps) 


xliQ 






XISO 
XISl 
XIS3 


0.8184 ± 0.0035 

1.556 ±0.005 
0.8327 ± 0.0035 


1.82 
3.09 
1.13 


1.543 ±0.029 
0.812 ±0.005 
1.490 ±0.011 


1.92 
0.69 
1.53 





was detected with Suzaku occurs at a phase of « 0.45. 
However, the discrepancy may be a reflection of longer 
term or stochastic variability effects unrelated to the 
rotation period of t Sco. 

It is also worth noting that all three detectors show 
rather high count rates at the last rotational phase. 
Deviations from the mean occur at the level of 3.4cr 
for XISO, 2.2cr for XISl, and 2.0cr for XIS3. The 
joint probability for this situation based on each light 
curve having independently a normal distribution of 
data points is « 10~ 7 . Interestingly, Figures|4}j6]do not 
indicate that the HR values are particularly anomalous. 
With all three detectors showing high count rates by 
amounts of 3.7%, 3.3%, and 2.0%, respectively, with 
essentially nominal values of HR, we suggest that the 
increased count rate is intrinsic to the source and es- 
sentially "gray" in character. Such a variation might 
result from an increase of hot plasma density or vol- 
ume filling factor. 

4. Discussion 

The critical conclusions from previous X-ray stud- 
ies are as follows. The X-ray lines of t Sco are strik- 
ingly narrow. The X-ray spectrum is dominated by a 
soft component rather typical of other hot stars. There 
is an abnormally strong hard component present as 
well, and this hard component in excess of 10 MK has 
been seen by ROSAT, ASCA, Chandra, XMM-Newton, 
and Suzaku. Moreover, the implied EM for the X-ray 
emitting plasma is around 5 x 10 54 cm 4 . As pointed 
out by Cassinelli et al. (1994), this competes with the 
EM available in the entire wind. Of course, it is not 
possible that the entire wind be at 1-10 MK since UV 
P Cygni lines of lower ions are observed. In fact, the 
problem is somewhat exacerbated by the expectation 
that in the standard wind shock model of line-driven 
winds, shocked hot plasma only begins to be gener- 
ated at 1 . 1-1 .5 stellar radii into the wind, implying that 



even less wind EM is available from which to generate 
the observed X-rays. This problem of large EMs for B 
star X-rays was highlighted in Cohen et al. (1997). 

In terms of generating the hard X-ray emissions, the 
magnetic field of t Sco is naturally an obvious candi- 
date. Although the occurrence of magnetism among 
OB stars is not yet known, there are efforts to de- 
termine the incidence rate (e.g., Bagnulo et al. 2006; 
Hubrig etal. 2008; Schnerr etal. 2008; Wade et al. 
2009). The impact of magnetism on X-ray emissions 
is also unclear. ud-Doula & Owocki (2002) used MHD 
simulations to study the influence of magnetism for the 
dynamics of line-driven winds from early-type stars. 
Naturally, the field topology takes on an increasingly 
dominant role when the magnetic energy density be- 
gins to exceed the ram pressure of the wind flow. 
When the magnetic energy density exceeds the wind 
kinetic energy density the field strongly modifies the 
wind flow leading to substantial volumes of wind con- 
finement. This sort of behavior was explored with a 
semi-analytic model by Babel & Montmerle (1997b). 

The strong field limit has scored successes in un- 
derstanding stars like 9 l Ori C and cr Ori E where 
X-ray data and UV and optical emission lines can be 
accurately modeled within a framework of a magnet- 
ically dominated wind flow. Less clear is how mod- 
est field strengths impact observables. Although Hen- 
richs, Schnerr, & ten Kulve (2005) have reported on 
correlations between periodic UV P Cygni line vari- 
ability with the incidence of magnetism, Petit et al. 
(2009) suggest that the X-ray properties of early-type 
stars are not so obviously correlated with the presence 
of surface magnetism. 

Our expectation is that the many unusual features of 
t Sco are somehow related to the presence of its mag- 
netic field. Therefore, the starting point for any anal- 
ysis should begin with the field topology. We possess 
a model for the magnetic field topology as extrapo- 
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Fig. 3. — A comparison of the areal sensitivities of 
the XISO and XIS3 detectors with XIS1. Top shows 
the XIS 1 detector response in red, whereas blue is for 
XIS0/XIS3. A percent difference of these is plotted in 
the middle panel. At bottom is an illustration of sev- 
eral spectra from XIS 1 (higher count rate in black) and 
XIS3 (lower count rate in green). The vertical dashed 
line demarcates the soft and hard bands used in our 
analysis. 

lated from the surface magnetograms acquired through 
Zeeman-Doppler imaging (Donati & Collier-Cameron 
1997; Donati etal. 1997; Donati et al. 1999; Donati 
et al. 2003). The extrapolation is performed using the 
"Potential Field Source Surface" method of Altschuler 
& Newkirk (1969) that has been used to extrapolate the 
Sun's coronal field from solar magnetograms. A code 
is used that was originally developed by van Ballegooi- 
jen, Cartledge, & Priest (1998). Since the method has 
been described in Jardine, Collier-Cameron, & Donati 
(2002), only a brief outline of the method is provided 
here. 

The magnetic field B is written in terms of a flux 
function *P such that B = -V*F, so the condition that 
the field is potential (i.e., V x B — 0) is satisfied auto- 
matically. The requirement that the field is divergence- 
free then reduces to Laplace's equation V~*¥ = with 
solution in spherical co-ordinates (r, i?, <p) given by : 



Fig. 4. — A plot of the hardness ratio as defined in the 
text plotted against the total count rate the XISO data 
of t Sco. The hashed region is the rms error ellipse for 
these data. A line connects the data points in chrono- 
logical and phase sequence, with #1 and #6 marking 
the first and last of the pointings (see Tab.[TJ. 



N 1 

V = J]J] [am/ + fci m r- (/+1) ] PU0) e^, (5) 

1=1 m=-l 

where the associated Legendre functions are denoted 
by P\ m . The coefficients a\ m and b\ m are determined by 
imposing the inferred radial field at the surface from 
the Zeeman-Doppler maps and by assuming that at 
some height /? s above the surface (known as the source 
surface), the gas pressure overcomes the ability of the 
magnetic field to confine it. Thus at the source sur- 
face, the field lines open up to become purely radial, 
and hence Bg(R s ) = 0. 

An illustration of the viewable field topology at the 
rotational phases of our pointings is shown in Fig- 
ure 13 Along the top, from left to right, are phases 
0.17, 0.34, and 0.50; at bottom the phases are 0.67, 
0.83, and 0.97 from left. In these figures the blue 
curves correspond to open magnetic field lines, and 
the white ones are closed loops. There is a dominant 
torus-like structure of closed field lines around an axis 
that is nearly perpendicular to the rotation axis. And 
on one hemisphere, there is a "mushroom-like" bundle 
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Fig. 5. — As in Fig. g] but for XIS1. Note that data 
points are too closely spaced to clearly mark the first 
pointing. 

of closed fields. This bundle feature is best perceived 
off the right limb for the middle bottom panel at phase 
4> = 0.83. At (p = 0.50 and 0.67, the "mushroom" lies 
most nearly forefront of the star, and nearly the entire 
torus region can be viewed. 

In some regions the magnetic field forms closed 
loops capable of confining plasma at coronal temper- 
atures, and in other regions the field lines are open to 
allow wind flow. At any given rotation phase, differ- 
ent field geometries are within view of a distant ob- 
server, leading to the expectation of rotational mod- 
ulation in the X-ray emission. If the X-ray emission 
were to come mainly from closed loops, then maxi- 
mum emission should occur when the largest number 
of closed loops are within view. Thus Donati et al. 
(2006b) predicted two X-ray minima at around phases 
0.3 and 0.8. For their coronal model, the overall mag- 
nitude of the EM at any given temperature depends on 
two parameters: the extent of the corona (which de- 
termines the magnetic geometry) and the plasma pres- 
sure at the coronal base (which determines the plasma 
structure). As an example, for a corona at 2 x 10 7 K 
and closed loops that extended no further than 1 stel- 
lar radius above the surface, they concluded that the 
closed loops alone could provide an emission measure 
of a few times 10 54 cirT 3 with a 40% rotational modu- 
lation. 



Fig. 6.— As in Fig.H but for XIS3. 

Our Suzaku XIS pointings reveal a spectrum that is 
predominantly soft, with a notable hard component, as 
has been consistently observed by other instruments. 
For the phases of 0.50 and 0.67, the "mushroom" of 
closed magnetic loops would be most nearly forefront, 
and yet somewhat surprisingly, no increase in hard X- 
rays was detected. Instead, we find that the XIS0 and 
XIS3 show little evidence for variability, but that the 
XIS 1 detector with about twice the sensitivity of the 
other detectors has a drop in total count rate near phase 
<p tss 0.5, with a 3.7cr deviation from the mean. From 
a comparison between XIS 1 spectra at different point- 
ings, there appears to be a slight deficit in soft emis- 
sion at <p « 0.5 at energies below 1 keV as compared 
to other pointings. This is suggested in Figure [5] where 
the lowest count rate for XIS 1 also has the highest HR 
ratio. We also note that all three detectors show an 
essentially "gray" increase in count rate at phase 0.98 
by about 3%, which corresponds to a time when the 
"mushroom" bundle of closed field lines is occulted 
by the star. 

Why is the hard emission not more variable? The 
HR values for XIS 1 are amazingly stable, and although 
HR values show greater dispersion for XIS0 and XIS3 
detectors than for XIS1, their variability is not sta- 
tistically especially significant. Our expectation was 
that the hard emission would form in the vicinity of 
the photospheric magnetic field, which asymmetrically 
distributed about this slowly rotating star. We there- 
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fore predicted that a rotational modulation could arise 
from stellar occultation. If the bulk of the hard emis- 
sion were to form at a characteristic radius Hmrd, then 
the geometric dilution factor 

W{f) = ^{\- (6) 

represents the areal fraction of a spherical shell of that 
radius that is occulted by the star, thus AQ occ /4^ = 
W(rh al -d)- If the hottest plasma forms characteristically 
at '"hard ~ 2/?* as indicated by the Chandra analysis of 
Cohen et al. (2003), then a fraction W ~ 7% of the 
source region would be occulted. 

If all of the X-rays were associated with the mag- 
netic field of r Sco and generated at 2R*, then given the 
distributed nature of the field reconstruction as shown 
in Figure|7] one might expect W to represent a charac- 
teristic scale of X-ray variability on the rotation phase. 
But this is a vastly oversimplified picture. In reality, 
t Sco has a complex field topology that is relatively 
strong which can alternately confine wind plasma in 
closed loops and modify wind flow in sectors of open 
field lines. There are a variety of emissive components 
at varying heights above the photosphere as a function 
of latitude and azimuth about the star. It is the recog- 
nition of these complexities that motivate the detailed 
field reconstruction presented in Jardine et al. (2002) 
and that was briefly summarized above. 

With its higher sensitivity, we interpret the mod- 
est variability in the total count rate of XIS 1 with rel- 
atively little change in the hard spectrum as a mo- 
ment constraint on the complex integral relations in- 
volving geometry and emissive components that de- 
termine the observed emission from an unresolved 
source. The overall relative absence of variability 
with phase should provide useful constraints for future 
modeling efforts that take into account both the wind 
flow on open field lines and confined plasma in closed 
field loops. 

In addition, a recent development is the identifica- 
tion of t Sco analogues. For a long time, the combina- 
tion of unusual properties displayed by r Sco had been 
unique. This may no longer be the case. The early 
B stars HD 66665 and HD 63425 had been noted to 
show unusual UV lines similar to r Sco (Massa, priv 
comm). Recently, magnetic fields have been detected 
in these two stars by the MiMeS collaboration^ (Petit 



2 www.physics.queensu.ca/ wade/mimes 



et al., in prep). It would be useful to learn if these two 
stars are also hard X-ray emitters. If so, then t Sco 
would no longer be a class unto itself, and the key to 
understanding r Sco may derive from studies of this 
set of three stars. 

5. Conclusions 

In Donati et al.'s (2006b) discovery paper of t Sco's 
magnetism, the authors predicted X-ray modulations 
on the rotation period of the star could be as much as 
~ 40%. In particular, two main eclipses were expected 
at phases of around 0.3 and 0.8 in conjunction with 
observed enhanced UV line absorptions at those times. 
Our pointings include phases close to those targets, but 
no variations anywhere near that amplitude were de- 
tected. It is possible that strong variability could have 
been missed if it had occurred in a relatively short in- 
terval of phase. However, rapid and high amplitude 
variability would appear to run counter to expectations 
based on the UV data. Donati et al. noted that an ab- 
sence of the predicted high amplitude X-ray modula- 
tion could be an indicator of smaller scale magnetic 
loops and confined hot gas across the stellar surface. 
Such loops would have evaded detection in their study. 

The most unexpected result of our study is that there 
is greater variability of the soft X-ray component than 
of the hard one. We nominally expected the soft emis- 
sion to be associated with wind shocks similar to OB 
stars, and the unusually bright hard emission to be as- 
sociated with the star's magnetic field. With closed 
loops extending out to a few stellar radii, variability of 
the hard emission was supposed to occur as a result of 
stellar occultation of the rotating magnetosphere; in- 
stead, the more sensitive XIS 1 detector shows a statis- 
tically significant low count rate event because the soft 
component drops more than the hard one. This came 
as a surprise. 

It thus appears that more effort will be needed to 
clarify the nature of t Sco's hard emission. Although 
the observed soft component is relatively nominal for 
a massive star wind of r Sco's early B spectral type, 
the standard wind shock theory most certainly does 
not predict large emission measures well in excess 
of 10 MK (e.g., Lucy 1982; Owocki, Castor, & Ry- 
bicki 1988; Feldmeier, Puis, & Pauldrach 1997). Our 
Suzaku study has not, in itself, resulted in a definitive 
explanation of how (or where) t Sco's hard emission is 
produced. There are a number of possible models that 
could be employed to interpret the new Suzaku data in 
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conjunction with data from other wave bands, such as 
the unusual UV P Cygni line shapes. We suggest that 
additional multiwavelength monitoring will be crucial 
for formulating a full picture of the circumstellar envi- 
ronment of r Sco. 

In particular, the phase coverage of our Suzaku 
study is still sparse. Plus it would be valuable to have 
multi-epoch data to determine whether the high count 
rates near phase 0.98 persists. In addition, the presence 
or absence of linear polarization would help to define 
the wind structure and/or the photospheric brightness 
distribution at inner wind radii. At IR wavelengths 
Brackett lines have been observed to be in emission 
whereas Her is in absorption for r Sco (Waters et al. 
1993). Although originally interpreted as a low density 
disk, Zaal et al. (1999) have suggested that the Brack- 
ett line emission could result from NLTE effects in the 
atmosphere of t Sco. A study of variability in these 
lines would help to determine their origin and place 
additional constraints on the density and temperature 



stratification that could be used in conjunction with the 
Zeeman maps and the X-ray and UV data to produce a 
more coherent picture of the r Sco environment. 
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